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Abstract 

Fluid-based locomotion at low Reynolds number is subject to the constraints of the scallop theorem, which 
dictate that body kinematics identical under a time- reversal symmetry (in particular, those with a single 
degree of freedom) cannot display locomotion on average. The implications of the theorem naturally compel 
one to ask whether similar symmetry constraints exist for locomotion in different environments. In this work 
we consider locomotion along a surface where forces are described by isotropic Coulomb friction. To address 
whether motions with a single degree of freedom can lead to transport, we analyze a model system consisting 
of two bodies whose separation distance undergoes periodic time variations. The behavior of the two-body 
system is entirely determined by the kinematic specification of their separation, the friction forces, and the 
mass of each body. We show that the constraints of the scallop theorem can be escaped in frictional media 
if two asymmetry conditions are met at the same time: the frictional forces of each body against the surface 
must be distinct and the time- variation of the body-body separation must vary asymmetrically in time (so 
quick-slow or slow-quick in the extension-contraction phases). Our results are demonstrated numerically 
and interpreted using asymptotic expansions. 



1. Introduction 

The capacity for locomotion is essential for the survival of much of life on Earth and is manifested in 
strategies as diverse as the organisms which depend on it. In fluids such as water or air and on scales 
ranging from microns to tens of meters , creatu res swim and fly by beating flagella, tails, wings, undulating 



their bodies, or actuating pumps (jVogell . 119941) . On land, animals crawl, walk, run, hop, climb, and slither 



using friction betw een their bodies and the ground, or undulate and appear to "swim" through sand or soil 



( Alexandej 



20031) 



The physics and scale of an environment determine the scope of successful propulsive strategies. The 
strategy employed by a scallop, for example, which is to quickly open its shell, displacing a large amount 
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of fluid, and then close it slowly, displacing a small amount of fluid, is ineffective if attempted in a fluid 
environment where inertial forces are overwhelmed by viscous forces. This notion forms the basis for the 
"scallop theorem" , which holds that locomotion at low Reynolds number is not possible if the kinematics of 
the body are identical under a time-reversal symmetry - which is always true if the deformation is controlled 
by a single degree of freedom (jPurcelll 119771) . The theorem relies on the linearity and time- independence of 
the equations of motion for the fluid (the Stokes equation) and states that in order to achieve self-propulsion, 
a body at low Reyn olds number must deform in a manner indicating a clear direction of time, for example 
in a waving motion (jLauga and Powersl 120091 ) . 



The beauty and important implications of this theorem, which greatly restrict the viable propulsive 
strategies available to organisms and machines at small Reynolds numbers, naturally compel one to seek 
similar fundamental results for other types of locomotion. For example, consider terrestrial locomotion 
for which organisms make use of friction forces between their bodies and a surface. Are strategies in this 
physical environment similarly constrained? 

In general, strategies for terrestrial locomotion may be divided into two categories: those relying on the 
movement of limbs (walking, hopping, climbing, and running) and those relying on movements of the body 
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2009). On crawling- types of locomotion, 



which is perhaps the simplest of all forms of terr estrial locomotion and is th e subject of this work, investiga 



tions have examined the locomotion of maggots (jBerrigan and Pepinl . 119951 ) and earthworms ( Quillin 



as well as hypothetical discre te-mass systems including s ystems consisting of two masses (Chernous'ko , 
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or kinematic constraints. 

A feature common to studies of the discrete-mass craw l ing strategies is that they consider system s 



with either many degrees of fre edom ( Zim merma n and et al 



anisotropic friction coefficients dZimmerman and et al 



2004, 



2007 



2009: 
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all 120041) , or analyze only locomotion in the case where 
the masses alternately stick and slide (jChernous Here we consider the most basic situation 

where deformation with one degree of freedom (and therefore time-reversible in its sequence of shapes) 
actuates a system with isotropic Coulomb friction coefficients and which always interacts with the surface 
through sliding friction forces, and propose to quantify the minimal requirements necessary to achieve 
locomotion. Is locomotion even possible with isotropic friction using time-reversible deformations? What 
are the associated minimum necessary mechanical or kinematic properties? How do they determine the 
direction and magnitude of motion? 




Figure 1: Two-body system (1 and 2) with prescribed relative distance, Z(t), translating along the x direction. 

To answer these questions, we consider a simple system consisting of two bodies which rest on a flat 
surface and are joined by a mechanism that enforces the time-variation of their separation distance. The 
variation of in the distance between the bodies gives rise to friction forces which determine their motion. 
Friction forces arc assumed to be isotropic (independent of the direction of body velocity) as it is evident 
that friction anisotropy will trivially lead to locomotion. Wc demonstrate that locomotion is possible in a 
system with only one degree of freedom provided that there is an asymmetry both in the friction properties 
of the system as well as the deformation kinematics (meaning the time-periodic variations in length display 
a quick-slow or slow-quick sequence of extension-contraction strokes). Finally, we develop a qualitative 
physical explanation for the mechanics o f strat egy of locomotion through an asymptotic analysis of the 



equations of motion (jBender and Orszaa 



1999T) and present a brief exploration of the space of physical 



parameters available to the system. 



2. Mathematical description of the two-body system 

Consider a system consisting of two bodies with position Xi (i = 1,2), velocity Xj, acceleration it, and 
mass rrii which are supported by a flat surface and arc constrained to translate along a horizontal line (the 
x direction). The bodies are connected by a mechanism which prescribes their relative distance and is the 
single degree of freedom available to the system to form the basis for its locomotion. The prescribed relative 
position of the bodies is denoted lit) = X2 — x\. A schematic representation of this system is shown in Fig.[TJ 

The variation in time of l(t) gives rise to a driving force, Fg,u which is exerted on each body by the 
linkage and is opposed by a friction force, Ff^, exerted on each body by the solid surface. The motion of 
the bodies is described by Newton's second law, 

1 



(F G ,i+F u ) 



(1) 



The friction force may be described by a variety of different models. Here, we restrict our analysis to 
isotropic Coulomb friction forces which are independent of the sliding velocity. We use a standard model 
which introduces a "static" friction force that requires that the force opposing it be greater in order for the 
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given body to translate, which we take to be the same as the sliding friction force. Mathematically we may 
express this by writing the friction law 

f -Fi a gn(xi) if lJb.il > F h 
{ ~F G ,i ii\F Gti \<Fi, 

where Fi is then the magnitude of the sliding friction force exerted on body i by the surface. Here we choose 
to represent the physical properties of the bodies and the surface as a friction force as opposed to a friction 
coefficient in order to simplify the analysis; friction coefficients fii can be found by computing /ij = Fijimg. 
Without loss of generality we assume that F\/m\ > -F2/7712 (or /ii > [12) to ensure that only body 1 sticks 
to the surface if the driving force is too weak to overcome the friction of both bodies. 

For the prescription of the relative distance between the bodies, we assume that l(t) is a periodic function 
which consists of an extension phase in which the bodies are pushed apart and a contraction phase in which 
the bodies are pulled together. Again without loss of generality, we specify the kinematics such that the 
extension phase is followed by the contraction phase and the bodies begin and end each period of the 
kinematic specification at their minimum separation. If the extension phase and contraction phase are 
opposite and equal in magnitude, we say that the kinematic specification is "symmetric" ; otherwise the 
kinematic specification is "asymmetric" . 

In order to determine the driving force between the bodies arising from the variation of their relative 
distance, we subtract the acceleration of each body to obtain 

x 1 -x 2 = l{t) = — (F f , 2 + F G . 2 ) - — {Ff ! + Fq 1) . (3) 
ni2 mi 

The relative position is enforced by a rigid mechanism which implies that Fq i = —F G ,2, and so we have 

F G2 = J^(l-^ + ^). (4) 
mi + m,2 \ m,2 mi J 

Inserting this relation into the equations of motion for each body, we find 

a, - 7712 u Fh + Ffl , 

m 1 +m 2 mi + m 2 ' 
mi F f2 +F h U 

mi + m 2 mi + m 2 

which completes our mathematical description of the body motion. For this equation of motion to hold 
for the blocks, the driving force must be large enough (and the kinematic acceleration large enough) such 
that both blocks slide. When the kinematic acceleration is too small for both blocks to slide, only the 
block with the smaller friction force will move. Also, since the kinematics are specified externally, the 
driving force is necessarily always large enough to force motion in one of the blocks. We therefore find 
that the Froude number, Fr, which is the ratio between driving forces and friction forces and is defined as 
Fr = (mi + m 2 )L/T 2 (Fi + F2), must be relatively large. In contrast, the locomotion of snakes is usually 
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associated with Fr < 1 (|Hu and et all 120091) , and while insufficient data is available to calcu late Fr for 



organisms which use closely related mechanisms (such as maggots, see lBerrigan and Pepinl (|l995l )). they are 
almost certainly associated with low Froude numbers as well. It seems likely that the reason for this is tied 
to the energy costs associated with maintaining high Froude numbers in frictional media. 

The most striking feature of the equations of motion is the discontinuous dependence of the friction force 
on the body velocity. This discontinuity is the critical nonlinearity which enables locomotion. However, 
while the equations are globally nonlinear, it is apparent that they are linear on intervals for which the 
direction of body motion does not change and that within these the position and velocity of the bodies can 
be found analytically. The full solution can therefore be determined by pasting together the exact solutions 
for adjacent intervals. Another technique to find steady solution computes the position and velocity of the 
system for a single period and equates initial and final conditions. 

3. Locomotion with one degree of freedom is possible 

With the mathematical description of body motion established, we ask whether sliding locomotion with 
a single degree of freedom is even possible. For our simple two-body system we observe two distinct modes 
of asymptotic motion, as illustrated in Fig. [2] stationary oscillation around a fixed point (no locomotion, see 
Fig. [2] A and B) and net translation (locomotion, see Fig.[5]C and D). In both examples from Fig. [21 bodies 
1 and 2 have the same mass (mi = m-i) but are subject to different friction forces on the surface (F2 7^ i*i). 
The two examples only differ in their kinematic specification: in A and B, the bodies are actuated with a 
symmetric kinematic specification, resulting in stationary oscillation, whereas in C and D an asymmetric 
kinematic specification leads to locomotion. In both cases the bodies exhibit transient behavior at the onset 
of motion which progresses into an asymptotic, stable mode of motion after long times. 

In general, three different regimes exist in the choice of relative friction force and kinematic specification: 
(a) equal friction forces and arbitrary kinematics, (b) unequal friction forces and symmetric kinematics, and 
(c) unequal friction forces and asymmetric kinematics. Regimes (a) and (b) lead to stationary oscillation, 
whereas parameter choices falling into regime (c) lead to translation. Locomotion with a single degree of 
freedom is therefore possible provided that there is an asymmetry in both the kinematics (actuation occurs 
at a different rate one way than the other) and the friction forces exerted by each block on the surface. 

4. Physics of locomotion with one degree of freedom: double symmetry-breaking 

The locomotion of the system depends on the frictional symmetry between the bodies as well as the 
symmetry between the extension and contraction phases of the kinematic specification. Physically, the 
asymmetry in friction forces is required because if the friction forces are symmetric, a solution exists in which 
the velocities of the bodies always oppose one another, and accordingly their friction forces cancel each other 
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Figure 2: Example of asymptotic stationary oscillation for equal masses (mi = ni2 = 1) but unequal friction forces (Fi = 4, 
F% = 1) and a symmetric kinematic specification. (A) Positions of both bodies (body 1, thick blue line; body 2, dashed red 
line) and geometric center (dash-dotted black line). (B) Time- variation of the normalized separation distance between the two 
bodies, (x2 — xi)/L, showing symmetric kinematics. (C) and (D): same as in (A) and (B) but with asymmetric kinematics 
leading to locomotion. 

such that their resulting equations of motion are merely proportional to the kinematic specification. The 
asymmetry in the kinematic specification is required because if the extension and contraction phases are 
opposite and equal in magnitude, a solution exists for which the net motion of the bodies over each phase 
is equal and opposite, which cumulatively result in no net motion over the kinematic period. The breaking 
of each of these symmetries divides the space of possible system properties into the three aforementioned 
fundamental regimes; we now examine each regime in detail. 

4-1- Symmetric friction forces 

When the friction forces exerted by each bodies are identical, the system does not translate. If the bodies 
are pushed apart in opposing directions after an initial time t = t$, the velocity and position of body 1 are 



G 



given by 

xi = ■ 1, (6) 

mi + m 2 

xi = ^r—[Kt)-l(to)]+xi(to), (7) 

mi + m 2 

where we have used ii(0) = 1(0) = 0. The position of body 2 can be found similarly, and if for simplicity 
we use xc(0) = 0, we find the position of the center of mass to be 

xc(t) = U^-l)[xi(t)-xi(t )}. (8) 

Since xi(t) is periodic in time (via Eq. [6]), we find that the position of the geometric center is periodic as 
well and the system exhibits stationary oscillation. When mi = m 2 , the center of the system is stationary 
for all times. 

4-2. Asymmetric friction forces, symmetric kinematics 

Further, even if the friction forces exerted by the bodies are asymmetric, the system will not translate if 
the kinematic specification is symmetric. This is best described mathematically by stating that l(t) = l(T—t) 
for all t, where T is the period of oscillation of lit). In this case, we are able to find a solution in which the 
system does not exhibit locomotion. To simplify the argument we consider the extension and contraction 
phases separately, which we denote A and B respectively. Each phase begins at a time to and ends at time 
tf, and the position of the geometric center during each phase is xa and xb- 

When the friction forces exerted by the bodies are unequal and the system is oriented as in Fig. [TJ we 
observe that the system velocity increases over the course of an extension phase. It is then possible to find 
a small negative initial system velocity at the beginning of the extension phase such that — iiyt(io) = &A{tf)- 
Then, by examination of the equations of motion ([TJ , it is evident that if the relative acceleration imposed 
by the contraction phase is opposite that imposed by the extension phase, then Xsito) = &A{tf) implies 
XB(to) = ~ x B(tf) and xa(^o) = x B(tf), and the considered initial system velocity constitutes an equilibrium 
initial system velocity in a periodic solution for the system motion. Further, the same logic applies to the 
system position: if the kinematic specification is opposite in sign and equal in magnitude, and if the body 
velocities are opposite in sign, then the motion of the geometric center over the extension and contraction 
phases, respectively, will be opposite in sign and equal, and the system will experience no net motion 
cumulatively over both phases. Locomotion is therefore not possible when the extension and contraction 
phases of the kinematic specification are symmetric, even under asymmetric friction. The statements of i j4.ll 
and t j4.2l constitute thus a form of scallop theorem for locomotion under isotropic friction. 

4-3. Asymmetric friction forces, asymmetric kinematics 

We have demonstrated that asymmetric friction forces and asymmetric kinematics can result in loco- 
motion, but the question remains: why? Some insight can be gained by considering the system in the 
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Figure 3: Illustration of body position (A) and velocity (B) in an extension phase when both bodies start from rest; (A): 
normalized body positions, Xj/L; (B): normalized body velocities, Xi X (T/L). There are two intervals distinguished by the 
signs of the body velocities; in the first, the bodies have opposite velocity and are pushed apart; in the second, both bodies 
translate forward. 

asymptotic limit of zero kinematic period, or fast kinematics and strong driving forces. We make two key 
observations. First, over the duration of an extension phase, the body with a smaller friction force (body 
2) will move farther forward (x > 0) than the body with a smaller friction force (body 1) when the bodies 
arc pushed part. This implies trivially that the geometric center of the system will move forward during an 
extension phase and backwards during a contraction phase. Second, if the extension phase is longer (and the 
kinematic velocity faster) than the contraction phase, the net change in system position over both intervals 
will be positive. This can be understood by examining the first terms of an asymptotic expansion in the 
extension phase duration for the change in system velocity and position over the phase. Due to the piecewise 
nature of the equations, we must consider separately the case where the bodies begin the extension phase 
at rest from the case where they begin with some initial velocity. 

4-3.1. A single extension phase when bodies are initially at rest 

Consider the two-block system when it is initially at rest at time to = 0. In this case, the extension 
phase is divided into two intervals distinguished by the sign of the velocity of body 1, as illustrated in Fig. [3] 
where we show both the body positions (A) and velocities (B). In the first interval, which constitutes the 
greater part of the extension phase, the bodies are pushed apart from their stationary position so that body 
2 moves forward and body 1 backwards, and l(t) > 0. This interval ends when the velocity of body 1, which 
is lesser in magnitude than body 2 due to its larger friction force, reaches zero and changes direction. In the 
second interval, both bodies move forward until the end of the interval is reached. 

The motion of each body is described by Eq. J5]). In the first interval we have that Ff 1 — Fx and 
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F 



h 



-F 2 , which yields 



Xl 



m 2 



-w 



F-i — F-2 



mi +m 2 mi + m 2 

Integrating from t = to t, and noting that the body is initially at rest, we find the body velocity, 

, \ m 2 ;, . Fi — F 2 

= r— *(*) + , *■ 

mi + m 2 mi + m 2 

and integrating again and setting Xi(0) = 0, we obtain the body position, 



xi(t) = — 



-l(t) 



(9) 



(10) 



(11) 



mi + 7Ti2 v ' m.i + m 2 2 

The first interval ends when the velocity of body 1 reaches zero. We expect this to occur very close to the 
end of the interval at a time t — T — AT, where AT <C T (see Fig. [3]). We therefore expand l(t) near t = T 
to find 



x\ (T - AT) = = AT m 2 l(T) - F x + F 2 + (Fi - F 2 )T + 0(AT 2 ) 



which implies 



AT = 



Fi — F 2 



-T. 



(12) 



(13) 



m 2 l(T) -F 1 +F 2 

To emphasize the dependence of quantities on the phase duration T, we note that if the length of extension 
remains constant we expect the kinematic acceleration to scale as I ~ 1/T 2 . We therefore introduce the 
notation l(T) = l T /T 2 , where l T is some constant which is proportional to the maximum separation between 
the bodies l(T), and write 

F, _ f„ _ 

(14) 



- 1 , 2 T 3 + 0(T 5 ). 



m 2 h 



The position of body 1 at the end of the interval is 



xi (T — AT) = — 



m 2 



-l(T) 



Fi - F 2 T 2 



0(T 3 



(15) 



mi + m 2 ' ' mi + m 2 2 

We can now find the final body velocities and positions. In the second interval we have that both body 
velocities are positive implying that Ff 1 = —F\ and Ff 2 = —F 2 . This yields the equation of motion 

m 2 »,.\ Fi + F 2 



-l{t) 



(16) 



mi+m 2 mi + m 2 

The velocity of body 1 is found by integration, where by specification the initial velocity of the body in this 
interval is zero, 



±i(t) = 

and at time ( = Twc find 



m 2 



mi + m 2 



l(t) - l(T - AT) 
Fi — F 2 



Fi + F 2 
mi + m 2 



[t — (T — AT)} 



±i(T) 



mi + m 2 



T + 0{T A ). 



(17) 
(18) 
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It is not necessary to solve for the velocity of body 2, because by definition of the kinematics it will equal 
the velocity of body 1 at time t = T . The position of body 1 is found by integrating the body velocity, 

Xl (T) = ^—l(T) + Fl ~ F2 + 0(T 4 ). (19) 

mi + m 2 mi + m 2 2 

Since x 2 (T) = xi{T) + l(T), we find that 

X2 (T) = mi l(T) +l Fl ~ F2 T 2 + 0(T% (20) 
y ' m 1 +m 2 2m 1 +m 2 

and the geometric center xc = (%i + x 2 )/2 is 

x c {T) = \^-^l(T) + J^-^T 2 + 0(T 4 ). (21) 
2 mi + m 2 2 mi + m 2 

The first term in this expression is the change in position due to the difference in body inertia and is 
independent of the phase duration. The second term, however, is proportional to the difference in friction 
forces between the blocks and therefore contains the key to the direction of body translation. Its origin is in 
the short interval of time that both bodies are moving forward: because an increase in phase duration causes 
the the length of this interval to increase, the change in final system position is correspondingly increased 
as well. 

4-3.2. A single extension phase when bodies have an initial positive velocity 

The situation detailed in the previous section is modified somewhat by the existence of an initial system 
velocity at the beginning of the kinematic period. A similar asymptotic analysis (the details of which are 
given in Appendix A) with an initial positive velocity xq yields the change in system velocity and position, 
respectively, 

xc(T) - ± = ^-^T - i ^i ( I -f)T* + 0(T 3 ), (22) 
mi + m 2 m 2 V'o h J 



and 



XciT) =2WT^- 2 liT) + X ° T+ 2WT^- 2 T ~ XO rnlT Q T + ° (T) ' (23) 



where It and Iq are lengths which characterize the acceleration term in the kinematic specification. The 
effect of the small initial velocity is to increase the total change in system position and to decrease the 
change in system velocity. The first is intuitive and easy to understand. The second is also intuitive and 
explains why the system exhibits stable dynamics, since period-wise increases in the initial velocity decay 
until the initial velocity reaches some steady-state. 

4-3.3. A full kinematic period when bodies translate in steady-state 

We can then address an extension phase followed by a contraction phase when the bodies are translating 
in steady-state. This requires us to first obtain the steady-state velocity for particular system parameters, 
and then insert this into our equation for the final body position after the full kinematic period. To derive 
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the steady-state velocity, we consider the change in velocity after an extension phase and a contraction phase 
and set the final and initial velocities equal to one another. 

For the contraction phase, the final velocity and position of the body when starting with some small 
initial positive velocity are simply the opposite of that for the extension phase, 

±c(T) - so = -^^T + ± ^i (± J- ) T 2 + 0(T% (24) 

mi + m 2 mi \/ , c h,c J 

and 

a; c (T = ■ l{T)+x T-- ■ T + x — T J + O T , (25 

2 mi + m 2 2 TOl + m 2 "ll/o,c 

where we have used the subscript c to differentiate the kinematic constants in the contraction phase (c) from 
the extension phase (e). The details of this calculation are found in Appendix B. For a contraction phase 
which is functionally equivalent and opposite to the extension phase, we have lo. e = — lo,c, and lx, e = —W.c, 
even if the phases have different durations. We then find the steady-state velocity by combining an extension 
phase with a contraction phase and and setting the final velocity equal to the initial velocity. By writing 
the extension phase duration as aT and the contraction phase duration as (1 — a)T with kinematic period 
T, we obtain 

. = (Fi - F 2 ) (2a - 1) / mim 2 \ 

2T{l/l -l/l T )[a 2 Fi mi + (l-a) 2 F 2 m 2 ] \m 1 +m 2 J [) ' [ ' 

We can simplify this expression by defining a reduced system mass m* = mim 2 / (mi +m 2 ) and lumping the 

kinematic-dependent parameters into a single parameter with dimensions of length L^in = Wt/[2(^t — lo)]] 

if L is the difference between the greatest and least separation of the bodies, then L^m — L for a piecewisc 

quadratic form for the kinematic specification and ir 2 L/8 for a sinusoidal form. The relation becomes 

x = (Fi-F 2 )(2a-l)^-= ^ k ;" m * , 2p ^+0(T). (27) 

a 2 Fimi + (1 - a) 2 F 2 m 2 T 

Accordingly we find that the steady-state initial velocity is inversely proportional to the total phase duration 
(or proportional to the prescribed kinematic velocity) . Using this initial velocity to compute the final position 
of the bodies, we find that 

x c = (Fi F 2 ) (2a - 1) — ^ k ;" m * + 0(T 4 ). (28) 

a z rimi + (1 — ay b 2 m 2 

This expression encapsulates the main result of this paper and explicitly establishes the two conditions 
necessary for locomotion. The first condition is that a 7^ 1/2; or there must be an asymmetry between the 
extension and contraction phases of the kinematics. The second condition is that F± ^ F 2 ; or there must be 
an asymmetry between the friction force exerted by each body on the supporting surface. It is important 
to note that condition requires an asymmetry in the friction force specifically, regardless of how this force 
depends on material parameters (such as mass or friction coefficient). Furthermore, we see that the change 
in position of the bodies therefore goes to a constant as the kinematic period vanishes, and the effective 
velocity of the system scales as xc/T ~ 1/T, or with the velocity of the kinematics. 
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Figure 4: Change in system position (A) and velocity (B) over a single extension phase versus phase duration; (A): change in 
normalized system position, Axc/L; (B) change in normalized system velocity, Axe x T/L. The full calculation is plotted 
with solid lines while the asymptotic expansion presented in the text for small phase durations is shown with dashed lines. 
The parameters are l(T) = 10, F2 = 1, mi = m2 = 1 and the plots are shown for three different friction asymmetries 
F 1 /F 2 = {2, 10,40}. 

4-4- A physical description for all parameter regimes 

With an explanation in hand for when the driving force is large compared to the friction force, we now 
examine the relationship between phase interval and the change in system velocity and position for the full 
range of extension phase durations. In Fig. 0] we plot the change in system position (A) and velocity (B) 
against the length of the extension phase for all physically valid extension phase durations and a kinematic 
specification which is piecewise quadratic and continuous in position and velocity. 

The dependence on the change in system position on phase duration (Fig. 0J\) provides an indication 
of the direction of system translation. We see that this dependence is monotonic, which implies that the 
system will always translate forward when the extension phase is longer than the contraction phase, and 
vice versa. 

The dependence of system velocity on phase duration (Fig. 0|3) provides an indication of the nature of 
the system dynamics. Here we find that the change in system velocity is not monotonic with phase duration. 
When the phase duration is short (large driving force, high Froudc number) the change in system velocity 
increases as the phase duration increases. When the phase duration is long (smaller driving force, relatively 
smaller Froudc number) this dependence is reversed. 

The combination of these two effects implies that there are two distinct regimes of translation for sliding 
locomotion of this system. When the change in system velocity is increasing with increasing duration, 
period-wise initial velocity of the system will be positive when the bodies are translating in the forward 
direction. However when the system velocity is decreases with increasing phase duration, we find that the 
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Figure 5: (A): Dependence of the period-averaged velocity, (U), on the kinematic asymmetry, a, for F2 = 1 and three values 
of Fi; (B) Dependence of {[/} on the friction asymmetry F1/F2 for a = 3/4 and three values of F2. In both cases we have 
mi = m 2 = 1, l(T) = 10. 



period-wise initial velocity of the system is antagonistic to the overall direction of translation. Finally, the 
maximum period-wise displacement of the system and therefore maximum average system velocity for a 
given kinematic specification occurs at the boundary between these two regimes. This can be explained by 
observing that a zero initial velocity at the beginning of the extension phase implies that the initial velocity 
for the contraction phase is at a maximum, which minimizes the decrease in system position over this phase. 

4-5. Parameter studies 

Once it is known that both an asymmetry in the specified kinematics and an asymmetry in the critical 
friction forces are necessary for the system to achieve net locomotion, we may ask how the motion of the 
system is dependent on the extent of the asymmetries. One metric which characterizes the capacity of the 
system for locomotion is its velocity, and so to probe the effect of physical asymmetry on system behavior 
we define a "period-averaged velocity" as 

x c (t Q + T) - x c (t ) 



(U) = 



T 



(29) 



For the kinematic specification, we choose a simple piecewise quadratic form and define a parameter a 
to characterize the asymmetry of the kinematics. The length of the extension phase is then aT and the 
contraction phase is (1 — a)T, and a — 1/2 implies a symmetric kinematic specification. 

In Fig. [51 we plot the period-averaged velocity against both the kinematic asymmetry (left) and the 
asymmetry in friction forces (right). For the kinematic asymmetry, we observe that the period-averaged 
velocity of the system increases as the specified kinematics become more asymmetric. 
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For the frictional asymmetry, we observe that the period-averaged velocity increases to a maximum and 
then decays to zero as the asymmetry increases, which is expected since no translation occurs cither for 
F1/F2 — 1 or as F1/F2 — > 00 (when the friction force is too large for the kinematic specification to force 
the system into a sliding mode of locomotion). Further, as explained in the previous section, we observe 
the two distinct regimes of translation in which the period-wise initial system velocity is cither in the same 
direction as translation (small F1/F2) or in the opposite direction (large F1/F2). 

5. Discussion 

In this work, wc attempt to understand the behavior of the conceptually simplest crawler possible: a 
one-dimensional, one-degree of freedom system consisting of two mechanically connected point masses. We 
find that such a system is able to achieve time-averaged translation even when its frictional interactions 
with the surface are isotropic and are mediated only through sliding contact, in addition to the more evident 
cases of intermittent static contact and anisotropic friction. We use computations and physical reasoning to 
show which and in what manner symmetries must to be broken to obtain net locomotion. In doing so, we 
demonstrate that: 

• Friction-based locomotion with one degree of freedom is possible because of the non-linear dependence 
of the friction force on body velocity; 

• Two symmetries must be broken in order for the system to achieve locomotion: the two components 
of the body must exert a different friction force on the supporting surface, and the body kinematics 
must be asymmetric on a single period of actuation, 

• The physical mechanism of locomotion results from an interval within a single kinematic period in 
which the body with a smaller friction force (faster) carries the other body (slower) forward 

• For each chosen set of physical parameters and kinematic function the system always eventually 
achieves a steady-state of locomotion; 

• Two fundamental regimes of translation exist which correspond to the relative strength of the driving 
force to the friction force (or Froude number) : when the driving force is large (higher Froude numbers) 
the system velocity at the beginning of each period is in the direction of translation; when the driving 
force is smaller (relatively lower Froude numbers) this initial system velocity opposes the direction of 
translation; and the maximum velocity occurs on the boundaries of these two regimes. 

Further, we observe that locomotion in this system occurs necessarily at relatively high Froude number. 
While experimental data on the forces exerted by crawling organisms is difficult to obtain, it seems likely 
that most crawling organisms locomote low Froude numbers. This may have something to do with the 
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energetic cost of maintaining high velocities (and therefore high Froude numbers) in frictional media; in our 
system, much mechanical energy is wasted in applying sufficiently large forces to the blocks to ensure that 
they always slide. 

The established criterion for the asymmetries of a two mass, friction-based system is fundamental knowl- 
edge in the physics of terrestrial locomotion. Additionally we hope that this study will help guide the 
analysis and design of simple crawlers for use in exploration and medicine, as well aid in the development 
of intuition and understanding of more advanced and complex modes of terrestrial locomotion. 
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Appendix A. Analysis of a single extension phase when bodies start with some small initial 
velocity 

We examine here the effect some small initial positive velocity has on the final position and velocity of 
the system after a single extension phase. When both bodies begin with a positive velocity the equation of 
motion for body 1 is 



TO2 ■■ F\ + F 2 
x\ = lit) (Appendix A.l) 

7771 +m 2 777l + "72 



Integration yields the body velocity, 



Tii F I F 

xUt) = tit) -t + ±0, (Appendix A.2) 

777l +777 2 777 1 + 7772 



and the body position, 



/ \ 771*2 / \ -^1 ~t~ F 2 t / A i 4 \ 

xAt) = lit) — + xot, (Appendix A. 3) 

777l + ■ r TL 2 mi + 7772 2 

where we set the initial position of body 1 to x±(0) = 0. The first interval ends at time t = t\ when the 
velocity of body 1 decreases to 0. In general this can only be solved if lit) is specified, but to start let us 
assume that t\ is small. We can then expand l{t\) to obtain 

777 F I F 

xAti) = — tj(0) -h + ±o + 0(t?), (Appendix A.4) 

777l + m 2 777l + "72 



We then find that x\(ti) = when 



777l + 7772 

ti = xri (Appendix A. 5) 

m 2 l(0)+F 1 +F 2 
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When T — > and using the notation 2(0) = lo/T 2 we have that 

i (mi +m 2 ) 2 ( 1 \ 
ti — ; 1 



m 2 l 

±o(mi + m 2 ) 2 
m 2 l 



J i y2 -F1+-F2 

7712^0 



■ (mi +M r2 , 
m 2 / 



m 2 l 

2 F 1 + F 2 
m 2 l 



m 2 l 



2<() 



(Appendix A. 6) 



0(T 6 ), 



which confirms our assumption that ii is small when T — >• 0. The position of body 1 at the end of the first 
interval is 



m 2 t 



F\ + F 2 t\ 
mi + Tthi 2 v ' mi + rri2 2 
ii m 2 Z(0) +Fi +F 2 \ 



Xoti, 



hi x - 



m\ + m 2 



(Appendix A. 7) 



x h 



0(T 2 ). 



In the second interval, body 1 is moving backwards and body 2 is moving forwards. The equation of motion 
is 



m 2 Ym , Fi-F 2 

X\ = : l{t) + : 

mi + m 2 mi + m 2 



Integration yields the body velocity 



xx(t) = 



m 2 



i(t)-i(tx) 



mi + m 2 

and integrating again yields the body position 

m 2 



Fx — F 2 
mi + m 2 



(t-tx), 



(Appendix A. £ 



(Appendix A.9) 



xx(t) 



mi + m 2 
, 1 Fx - F 2 



[l(t)-l(tx)} 



m 2 



mi + m 2 



l(tx)(t-tx) 



{t-txf + xx{tx). 



(Appendix A. 10) 



2 mi + m 2 

This interval ends at time t = t 2 when the velocity of body 1 goes to zero again and both bodies are 
translating forwards. We use the ansatz that t 2 is very close to the end of the interval T and, as in the case 
where the bodies start from rest, define AT = T — t 2 . We then have that 



l(T - AT) = i(T) - ATZ(T) + 0(AT 2 ), 
= -AT[(T) + 0(AT 2 ), 



(Appendix A. 11) 



and 



l(T - AT) = 1{T) - AT/(T) 



AT 5 



l'(T) + 0(AT 3 



AT 5 



(Appendix A. 12) 



l(T) + 0{AT 3 ). 
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The equation for body velocity becomes 
m 2 



= 



mi + m 2 



-ATl(T) — *ii(0) 



F\ — F 2 
mi + m 2 



[(T - AT) - h) 



AT 



m 2 l(T) ~ {Fi - F 2 ) +h m 2 i(0) — (F± — F 2 ) 



+ (Fi - F 2 )T, 



AT 



m 2 l(T) - (Fi - F 2 ) + x a (mi + m 2 ) 



m 2 i(0)-Fi+F 2 
m 2 l(0) + Fi+F 2 



This implies that 
AT 



(Fi - F 2 )T, 



1 



m 2 l(T) -Fi+F 2 
■■■ + (Fi- F 2 )T 



x (mi + m 2 ) 



m 2 j(0) -Fi+F 2 
m 2 l(0) + Fi+F 2 



Using the notation [(0) = l /T 2 and l(T) = l T /T 2 , we find 
T 2 / 1 



AT 



i 2 l T \ i _ t2£l 



m 2 'o 



(Fi - F 2 )T 



T 



.1 + T' 

m 2 lq 

: x (mi + m 2 ) 
m 2 lx 
2 Fi+F 2 



m 2 lj 



1 + T 



0(T 7 ). 



Fi — F 2 
m 2 l T 



1 - T 



\F\ — F 2 
m 2 l 



1 -T 



m 2 l 



o 



0(T 6 ). 



Rearranging, we find 
AT = 



Xo(rni + m 2 ) 



t2 _ i Fi - F 2 i ^ 3 
m 2 lT 



m 2 l T 

io{mi + m 2 ) ( Fi - F 2 2f\ 



m 2 lx 



m 2 l T m 2 l 



T 4 + 0(T 5 ). 
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The position of body 1 at the end of the second interval is 

m 2 



x\{T — AT) = - 



mi + 7712 

m 2 



+ 



mi + m 2 

1 F i - F 2 

2 mi + 7712 



mi + m 2 
m 2 



p(T- AT) -/(*!)] 
i(h)(T - AT - h) 
(T-AT-ti) 2 +xi(h), 

AT 2 l(T) - t 2 J(0) 



mn , ,„ N 1 m? 
l(T) - 



+ 



mi + m 2 
1 Fi-F 2 



2 mi + m 2 
tJ(0)(T - AT - h) 

2 , x h 



2 mi + rrt2 
m 2 



(T — AT — ti] 



0(T% 



mi + ?ri2 



l(T)+T*[- 



1 Fi - F 2 



2 mi + rri2 



AT 



1 m 2 /(T) 1 Fi - T 2 



2 mi + m.2 2 mi + m2 

1 m 2 f(0) | 1 F - F 2 

2 mi + m-2 2 mi + m,2 



TAT 



hAT 



tiT 



Fi — F 2 
mi + m2 

m 2 /(0) [ Fi-F 2 
mi + m2 mi + m 2 



m 2 l(0) 



Fi — F 2 



mi + m 2 mi + m 2 



0(T 4 ). 



After much algebra, we obtain 

xi(T - AT) = - 



m 2 



mi + m 2 



l(T) + x T 



+ 



1 



Fl — F2 
mi + m 2 



(*o) 



2 mi + m 2 



m 2 lq 



xo_ f Fi-F 2 _ 2F\ 
m 2 \ l T l 



j,3 + Q(j,4 



(Appendix A. 17) 



(Appendix A. 18) 



In the final interval both bodies have a positive velocity. The equation of motion for body 1 is 

(Appendix A. 19) 



mo . Fi + F 2 
-l(t) - 



mi + m,2 mi + m 2 

so, noting that body 1 is at rest at the beginning of the final interval, the velocity of body 1 is 



mo 



mi + m 2 



l(t) - l(t 2 ) 



Fi + F 2 
mi + m 2 



(t-ta), 



(Appendix A.20) 



and its position, 



Xi 



m 2 



mi + m 2 

1 Fi+F 2 

2 mi + m 2 



[l(t] - l(t 2 )} 



m 2 



mi + m 2 



-i(h)(t-t 2 ) 



(t-t 2 ) 2 + Xi(t 2 ). 
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(Appendix A.21) 



At time t = T wc find for the body velocity 

m 2 



±i(T) 



mi + m 2 

m 2 



/(T) - i(T - AT) 



F 1 +F 2 
mi + m 2 



AT, 



-AT 
-AT 



mi + m 2 
m 2 



mi + m 2 



Z(T)- 
Z T T" 



Ti + T 2 
mi + m 2 

2 , F l +F 2 



mi + m 2 



0(T 4 ), 
+ 0(T 4 ), 



(Appendix A.22) 



. . Fl — F 2 . 2F\ / 1 1 \ t2 

x H ■ T - xo 7- - -r- T + G>(T^ 

mi + m 2 m 2 V to <t/ 



and for the position of body 1 
m 2 



£1 



[Z(T) - Z(T - AT)] 



m 2 



mi + m 2 

- - Fl+f2 AT 2 + xi(T - AT) , 
2 mi + m,2 

AT 2 Z(T) - - 



mi + m 2 



-Z(T - AT)AT 



2 mi + m-2 

1 m 2 Z T AT 2 

2 mi + m2 T 2 

m 2 



2 mi + m 2 
a?i(T- AT) + 0(T 4 ) 



AT 2 + xx{T - AT) + 0(T 6 ), 



(Appendix A.23) 



l(T)+x T- 



1 / Fi 



. 2Fi 



T 2 + x -T A + 0(T 4 ). 



mi + m 2 " ' 2 \ v mi+m 2 y m 2 io 

At the end of the extension phase, the velocity of body 2 is equal to that of body 1, and the position of 
body 2 is x 2 (T) = x\(T) + l(T). The velocity of the geometric center ±c = {x\ + i 2 )/2 is therefore 



i c {T) = x + I^J±t - xo^ (r-r) T2 + °( t3 )' 

mi + m 2 m 2 \l It ' 



and the position of the geometric center xc = {x\ + x 2 )/2 is 



xc(T) = - 



1 mi — m 2 



2 mi + m 2 



?(T) + i T+- 



1 / Fi - F 2 



2 V mi + m 2 



. 2Fi 



7 1 +±o ^T J + 0(T 4 ). 

m 2 i 



(Appendix A.24) 



(Appendix A.25) 



Appendix B. Analysis of a single contraction phase when bodies start with some small initial 
velocity 

We examine here the change in velocity of the two body system after a contraction phase when the 
bodies begin and end with some small initial positive velocity This velocity is presumed small enough such 
that the contraction kinematics induce a reversal in the velocity of body 2, but large enough such that the 
body still has a positive velocity at the end of the phase. This analysis is very similar to the extension phase 
analysis, except that it is the velocity of body 2 which is now of interest. The equation of motion for body 
2 is 



x 2 



mi 



mi + m 2 



i(t) 



Fl + F2 
mi + m 2 



(Appendix B.l) 
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This equation of motion holds until the velocity of body 2 changes sign. Using the same steps as for the 
extension phase and letting T — > 0, we find the time h at which x 2 (t\) = 0, 

mi + m 2 



h = Xq- 



mil 

The position of body 2 at t — h is 



iJ(Q) +F 1 +F 2 ' 
mi + m2 T 2 (l + T 2 ^±^)+0(T% 



(Appendix B.2) 



x 2 (h) = x 2 (o) + — (i(h) - i(o)) - I 1 2 1 2 + x h, 



x 2 (0)+t 



mi + 1^12 v v *' " " 2 mi + ?7i2 
*i Fx + F 2 - l'(0)mi 



a; 



mi + m 2 



(Appendix B.3) 



x 2 (0) + ^ + O(Ti). 



In the second interval, body 1 is moving forwards and body 2 is moving backwards, and the equation of 
motion for body 2 is 

x 2 = lit) -• (Appendix B.4) 

mi + m 2 mi + m 2 

As in the extension phase with body 1, this interval ends when the velocity of body 2 changes sign. It is 
convenient to again define AT = T — t 2l and we find 



AT 



mi l(T) -F 1 +F 2 
- (Ti - F 2 )T 



xo (mi + m 2 ) 



mj(0) ~Fi+F 2 
mj(0) -Fi-F 2 



(Appendix B.5) 



With the notation 1(0) = Iq/T 2 and l(T) = It/T 2 , and expanding in powers of T, this becomes 



AT 



Ti - F 2r ^ 3 f +T 2 F i- F 2 



mih 



0(T 7 ) 



T2 x (mi +m 2 ) f +T - 2 Fi- F 2 



1 + T 



milx 
2 Fi + F 2 



mxlo 



ttliIt 
+ 0(T 6 ), 
Fi-FA t3 

mi^T V m i^T / 

x (mi + m 2 ) ( Fi - F 2 2F 2 



1 - T 



2 Fi - F 2 
mil Q 



(Appendix B.6) 



x (mi + m 2 ) 



T 2 



ttiiIt 



milx mild 



T 4 + (9(T 5 ). 
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Expanding the kinematic specification in the same way as in the extension phase calculation, and noting 
that in this case l(T) = 0, we find the position of body 2 at the end of the second interval, 

'"l 



x 2 (T — AT) 



mi + m 2 
mi 



mi + m 2 
1 Fi - F 2 



[l(T- AT) - l(ti)] 
i(h) (T - AT - h) 
(T-AT-tif+x 2 (ti), 



2 771 1 + ni2 
TUi ,,„., 1 



mi + m 2 
nil 



1(0) 



in i 



mi + 7712 

1 Fi - F 2 



2 mi + 7712 

tj(0) (T-AT-ti 

2 



AT 2 l(T) - t 2 J(0) 



2 mi + 7712 
7«i 



(T- AT- h) 



1 Fi-F 2 



mi + m 2 



1(0) - f - 



tiT 



AT 2 



2 mi + m 2 
Fi — F 2 



mj(0) 
mi + m 2 mi + m 2 

1 777! l(T) 1 Fi - F 2 



x 2 (0) 



TAT 



i ti 



0(T 4 ), 



Fi — F 2 
mi + m 2 



2 777l + "72 2 777l + m 2 

mi 1(0) \ Fi-F 2 \ 



hAT 

x 2 (0) 



mi + 7772 2 mi + m 2 I 
mj(0) Fi - F 2 



mi + m 2 
x ti 



mi + 7772 

0(T% 



which simplifies to 



x 2 (T - AT) = x 2 (0) - 



nil 



mi + m 2 



1(0) + x Q T 



1 
2 

x 
mi 



Fl - F 2 2 7771 + 777 2 

+ (x ) 



mi + m 2 
Fi — F 2 
It 



2F\ 
lo 



T 3 + 0(T 4 )- 
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In the third interval, both bodies are translating forward. The equation of motion for body 2 is then 



x 2 



mi 



mi + 7772 



l(t) 



Fi + F 2 
mi + m 2 



(Appendix B.9) 
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By integrating from T — AT to T and employing the fact that x 2 (T — AT) = 0, we find the velocity of body 
2 at the end of the contraction phase, 



x 2 (T) 



mi 



mi + m 2 
nil 



l(T) - l(T - AT) 



F± + F2 



mi + m 2 



AT, 



mi + m 2 

Fi — F 2 



ATl{T) - Fl + F2 AT, 



x 



mi + m 2 



mi + m 2 
2F 2 



(Appendix B.10) 



T + x Q 



mi \l 



1 

h 



The position of body 2 is 



*2 (T) 



11} 1 



mi + m 2 

1 Fi+F 2 



[l(T) - l(T - AT)} 



mi 



mi + m 2 



0(T 3 ). 



4(T - AT) 



2 mi H 
mi 
2 mi + m 2 



m 2 



AT 2 + x 2 (T — AT), 



1 



AT 2 /(T) + x 2 (T - AT) + 0(T 4 ), 



(Appendix B.ll) 



x 2 (0)- 
+ x Q 



in 1 



mi + m 2 
2F 2 



1 Ft - Fo 
l(0) + x T---^— 

I mi + rri'2 



T 



mil 



T 3 + 0(T 4 ) 



Finally, we find the velocity of the geometric center by observing that at the end of the contraction phase we 
have ii(T) = x 2 (T) = xc(T), and we find the position of the geometric center by observing that because 
1{T) = then we must have xc{T) = x 2 (T). Further, if we presume that, as in the extension interval 
calculation, that the initial position of the geometric center is xc(0) = 0, we find that x 2 (0) = Z(0)/2. The 
velocity of the geometric center at the end of the contraction phase is therefore 



x c (T) = ±0 



( Fi — F 2 



\mi + m 2 

and the position of the geometric center is 



mi W It' 



xc{T) 



1 mi — m 2 



2 mi 



m 2 



l(Q) + x T 



Z mi + m 2 iTiik) 



(Appendix B.12) 



(Appendix B.13) 
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